Relaxation dynamics of charged gravitational collapse 
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We study analytically the relaxation dynamics of charged test fields left outside a newly born 
charged black hole. In particular, we obtain a simple analytic expression for the fundamental 
quasinormal resonances of near-extremal Reissner-Nordstrom black holes. The formula is expressed 
in terms of the black-hole physical parameters: u) = g$ — i2TTTBH{n+ |), where Tbh and $ are the 
temperature and electric potential of the black hole, and q is the charge of the field. 
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The radiative test fields of complete gravitational collapse decay with time leaving behind a 'bald' black hole. This 
is the essence of the no- hair conjecture introduced by Wheeler in the early 1970's The various no- hair theorems 
, state that the external field of a dynamically formed black hole (or a perturbed black hole) relaxes to a Kerr-Newman 
spacetime, characterized solely by the black-hole mass, charge, and angular momentum. This implies that test fields 
left outside the newly born black hole would either be radiated away to infinity or swallowed by the black hole. 
The relaxation phase in the dynamics of perturbed black holes is characterized by 'quasinormal ringing', damped 
, oscillations with a discrete spectrum (see e.g. Q for a review and detailed lists of references). At late times, all 
^sj ' perturbations (test fields) are radiated away in a manner reminiscent of the last pure dying tones of a ringing bell 
', [sHa]. The black hole quasinormal modes (QNMs) correspond to solutions of the perturbations equations with the 
^ ' physical boundary conditions of purely outgoing waves at spatial infinity and purely ingoing waves crossing the event 
horizon. Such boundary conditions single out a discrete set of black- hole resonances {uJn} (assuming a time dependence 
^ ■ of the form e"^'^*). 

{3jr)| In accord with the spirit of the no-hair conjecture the external test fields would either fall into the black hole 
or radiate to infinity. This implies that the perturbation (test field) decays with time and the corresponding QNM 
I frequencies are therefore complex. It turns out that there exist an infinite number of quasinormal modes, characterizing 
^ oscillations with decreasing relaxation times (increasing imaginary part), see 0^ and references therein. The mode 
0^ ' with the smallest imaginary part, known as the fundamental mode, determines the characteristic dynamical timescale 
m for generic perturbations to decay [lol - [l6| . [We recall that, in asymptotically flat spacetimes the quasinormal ringing 
is followed by an inverse power-law decaying tail. However, it should be noted that most of the energy of the test field 
is radiated away to infinity (or swallowed by the black hole) in the quasinormal ringing phase [2]. The decaying tail 
contains only a small fraction of the field's energy. It is therefore physically justified to determine the characteristic 
lifetime of the perturbation (of the test field) from the decay of its associated QNMs.] 

In this work we study the nearly spherical gravitational collapse of charged matter to form a charged Reissner- 
Nordstrom (RN) black hole. The black-hole formation is followed by a relaxation phase which describes the decay of 
two types of perturbation fields (test fields): (1) coupled gravitational-electromagnetic test fields, and (2) charged test 
fields. In the present work we focus on the relaxation dynamics of the charged test fields left outside the newly born 



' black hole. In particular, we shall determine the fundamental (least-damped) resonant frequencies of a charged scalar 
- - 1 field in the charged black-hole spacetime. As we shall show below, the spectrum of charged quasinormal resonances 
can be studied analytically in the near-extremal limit Q — >■ M (we use natural units for which G = c = h — 1), where 
M and Q are the mass and charge of the black hole, respectively. [This also corresponds to the Tbh ^ limit, where 

(M^ _q2\1/2 

Tbh = 2^[jJ^(jv/2_Q2-)i/2]2 is the Bekenstein-Hawking temperature of the black hole.] 

The dynamics of a charged, massive scalar field vj/ in the RN spacetime is governed by the Klein-Gordon equation 

[{\/''~iqAn{^,~iqA,)^fi^]^^0 , (1) 

where Ai, is the electromagnetic potential. Here q and /i are the charge and mass of the field, respectively, [q and fi 
stand for q/h and n/h respectively, and they have the dimensions of (length) ^^.] One may decompose the field as 

^i,n{t,r,e,4,) = e'"'^Si„,{e)Ri„^{r)e-'^' , (2) 

where {t, r, 6, (p) are the standard Schwarzschild coordinates, w is the (conserved) frequency of the mode, I is the 
spherical harmonic index, and m is the azimuthal harmonic index with < m < I. (We shall henceforth omit the 
indices I and m for brevity.) With the decomposition©, R and S obey radial and angular equations both of confluent 



Heun type coupled by a separation constant A [171 [1 
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The angular functions S{0) are the famihar spherical harmonics which are solutions of the angular equation 18 1 

sm0de\ 89 J V sin^6'/ 

The angular functions are required to be regular at the poles 6 = and = tt. These boundary conditions yield the 
well-known set of eigenvalues Ai = Z(Z + 1), where I is an integer. 
The radial Klein- Gordon equation is given by [l9| 



dr ( dr ) 



- A^i^r-" + R = , (4) 



where A = — 2Mr + Q"^ and K = ujr'^ — qQr. The zeroes of A, r± = M ± {AP — Q'^Y^, are the black hole (event 
and inner) horizons. 

In order to determine the free resonances of the field one should impose the physical boundary conditions of purely 
ingoing waves at the black-hole horizon (as measured by a comoving observer) and purely outgoing waves at spatial 
infinity. That is, 
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I g-»((^-<7<E.)y as r^r+ (y -> -cx)) 



where the "tortoise" radial coordinate y is defined by dy = (r^/A)(ir and $ = Q/r+ is the electric potential at the 
horizon. These boundary conditions single out a discrete set of resonances {w„} which correspond to free oscillations 
of the field. These resonances determine the ringdown response of a black hole to external perturbations. Note that 
the boundary conditions of outgoing waves at spatial infinity imply a; > /i. 
It is convenient to define new dimensionless variables 

_ r+ - r_ _ uj-q<P 

X = ; T = ; m = — — ; k = 2ujr+ - qQ , (6) 

r+ r+ 2ttTbh 

in terms of which the radial equation becomes 

x{x + t)^ + {2x + t)^ + VR = , (7) 
dx'^ dx 

where V = K"^ /r\x{x + t) — [/i^r5|_(x + 1)^ + l{l + 1)] and K = r']^Lux'^ + r^kx + r^vDT/2. As we shall now show, the 
radial equation is amenable to an analytic treatment in the regime r ^ 1 with M{lj — (7$) <C 1. 

We first consider the radial equation ([7]) in the far region x 3> maxjr, Af (a;— (7$)}. Then Eq. ([T]) is well approximated 

by 

^'0 + 2x^ + Ff,ri? = O, (8) 

where Vfar = (w^ — /i^)r^a;^ + 2{Lok — fPr+)r+x — [fi'^r\ + l{l + 1) — fc^]. A solution of Eq. ([8]) that satisfies the 
boundary condition ([5]) can be expressed in terms of the confiuent hypergeometric functions M{a, b, z) 2^ 21 1 

R ^ (7i(2iV'Jj2— ^^^)i+.52.-H«e-'V^^^^'^+^M(i +i5 + iK,l + 2i5, 2i^J^^^r+x) + C2(5 ^ -8) , (9) 



where Ci and C2 are constants and 

S^^k'- M^r^ - (/ + 1)2 ; (^0) 



The notation {S — —S) means "replace 5 by —5 in the preceding term." 

We next consider the near horizon region x 1. The radial equation is given by Eq. ([7]) with V — >■ Vncar = 
— [^^r^ + l{l + 1)] + (kx + wt/2)'^/x{x + r). The physical solution obeying the ingoing boundary conditions at the 



horizon is given by [20|, |21 1 



1)] + (k 
■by 0, 



i fx \^{\^^^) 1 1 

i? = a;-2^(_ + 1) ^FY{^^i6-ik,-~i8-ik;\-iw;-xlT), (11) 

\t / 2 2 
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where 2^1(0, b; c; z) is the hypergeometric function. 

The solution s (pi) a nd pT|) can be matched in the overlap region niax{r, M{uj — q^)} ^ a; ^ 1. The x ^ 1 limit of 
Eq. dH) yields [IJIII 



The a; > r limit of Eq. (HH) yields [13, 1!! 



r(2i(5)r(i - iro) 



r(i + i(5 - ifc)r(i + j(5 - iro + z/c) 
By matching the two solutions in the overlap region one finds 



-iS—itu/2 



r(2i(5)r(l - im) 



r(i + iS- ik)T{^ +iS-iu7 + ik) 



— ^i^r^) 



r(-2i(5)r(i - iro) 



— i5 — ik)T(=2 ~ iS — iw + ik) 
Approximating Eq. © for a; 00 one gets [2^, [2l| 



-i(5 



(12) 

(13) 

(14) 
(15) 



R 



Ci(2zv/^J^r+)^^-^ii±^x-i+^'= + C2(<5 ^ -5) 
1 (5 + id + zk) 

^ ,o- — 1 r(l + 2z5) 

C,(2.V^^r,) ^^r^-^^- 



(16) 



A free oscillation (a quasinormal resonance) is characterized by a purely outgoing wave at spatial infinity. The 

coefficient of the ingoing wave e^*v^^^''+^ in Eq. (jl6p should therefore vanish. Taking cognizance of Eqs. (jl4p - 
([TC]) . one finds the resonance condition for the free oscillations of the field: 



r(2i(5)r(l + 2iS){-2iTy/uj^ ~ fi^r+y 



r(-2i(5)r(l - 2i5){-2iTy/uj^ - ii^r+) 



+ iS — iK)r(i + iS — ifc)r(i + i6 — iw + ik) ^ — zK)r(i — i6 — ik)T{^ — iS — iw + ik) 



= 



(17) 



The resonance condition (IT71) can be solved analytically in the regime r ^ 1 with w ~ g<i>; We first write it in the 
form 



1 



r(i — i(5 — inj + ik) 



V X (-2iTv/a;2 - ^2^+) 



-2i5 



(18) 



where V = [T{2i5)fT{\ - iS - m)r(i - iS - ik) / [T {-2id)]^T +iS- iK)T{^ +iS- ik)T{^ + iS-iw + ik). We note 
that T> has a well defined limit as Q — ?> M and cj — > (7$. 

In the limit w — > (7$, where to is almost purely real, one finds from Eq. (flUl) that (52 is also almost purely real. 
As we shall now show, the resonance condition Eq. (fT8|) can be solved analytically if (5 ^ 1. Thus, we restrict our 
attention to the regime > ^? + 0{M''^). Then one has {-i)''^'-^ = e(~*5)("245) = g-^r* ^ Qne therefore finds 
e = (-2iTv'a;2 - ^2j._^)-2«5 ^ 

Thus, a consistent solution of the resonance condition (fT8|) may be obtained if l/r(i — z(5 — iro + i/c) — 0(e) [2^ . 
Suppose 



— — i(5 — 1-07 + ik = —n + r/e + O(e^) , 
where n > is a non-negative integer and 77 is a constant to be determined below. Then one has 



r(— — iS — iw + ik) ~ r(— n + rje) ~ {~n) ^r{—n + 1 + ?/e) 



[{~l)-nir'r{rje) 



(19) 



(20) 



where we have used the relation r(z + 1) = zT{z) [2^. Next, using the series expansion l/r(z) — J^kLiCkz'^ with 
ci = 1 [see Eq. (6.1.34) of [13]], one obtains 



l/r(2 -iS-iw + ik) = (-l)"n!77e + 0(6^) 



(21) 
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Substituting (UJ) into ^ one finds 77 = 2?/[(-l)"-n!]. 

Finally, substituting w = {oj - q^)/2nTBH and k = qQ + 0{MTbh) for w = (7$ + 0{Tbh) into Eq. (HH), one 
obtains a simple formula for the quasinormal resonances of the charged field: 

^ = — - i27:TBH[{n + ^ - V^) - - qQ)] , (22) 
r+ 2 

where <5 = - ^2)q2 _ _^ l/2)2]i/2 + 0{MTbh) and n = 0, 1,2, ... . 

As emphasized above, our boundary conditions require > /i. From Eq. ([22l) one finds SRcj = (7[1 + 0{MTbh)\ 
in the near-extremal limit MTbh — >■ 0. Thus, our analytical treatment is valid in the regime qj yi.> 1 + 0{MTbh)- 

Taking cognizance of Eq. (122]) one finds Qlu = 0{Tbh) for the charged test fields in the near-extremal limit. We note 
that numerical studies have indicated that neutral test fields (coupled gravitational-electromagnetic perturbations and 
Dirac perturbations) are characterized by 9w = 0{M~^) ^j. One therefore realizes that the relaxation dynamics of 
the newly born charged black hole is dominated by the charged perturbation (test) fields (Note that Tbh ^ M^^ in 
the near-extremal limit.) We would like to emphasize that our analytical approximation is valid for real values of 5 
with 5^1. This implies that we cannot extrapolate our results to the case of neutral test fields: Had we taken in 
Eq. (jlOp 3?w ~ (7$ — >■ for the neutral test fields, we would have found that i5 is imaginary in this limit- this would 
be beyond the regime of validity of our analytical approximation. 

In summary, we have analyzed the relaxation dynamics of charged gravitational collapse. In particular, we have 
studied analytically the quasinormal mode spectrum of charged fields in the spacetime of a near-extremal charged 
black hole. It was shown that the fundamental resonances can be expressed in terms of the black-hole physical 
parameters: the temperature Tbh and the electric potential $. 

It is worth mentioning that a fundamental bound on the relaxation time Tioiax of a perturbed thcrmodynamical 



system has recently been suggested [llMld |24|. Treiax > Tl/ttT, where T is the system's temperature. This bound can 
be regarded as a quantitative formulation of the third law of thermodynamics. Taking cognizance of this relaxation 
bound, one may deduce an upper bound on the black-hole fundamental (slowest damped) frequency 

min{3a;} < ttTbh/^ ■ (23) 

Thus the relaxation bound implies that a black hole should have (at least) one quasinormal resonance whose imaginary 
part conform to the upper bound (|23|1 . This mode would dominate the relaxation dynamics of the newly born black 
hole and will determine its characteristic relaxation timescale. Taking cognizance of Eq. (f22|) and substituting 
n = for the fundamental quasinormal resonance, one obtains minj^w} — ■kTbh/^- One therefore concludes that 
near-extremal charged black holes actually saturate the universal relaxation bound. 
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